Jamming of Soft Particles: Geometry, Mechanics, Scaling and Isostaticity 
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Amorphous materials as diverse as foams, emulsions, colloidal suspensions and granular media 
can jam into a rigid, disordered state where they withstand finite shear stresses before yielding. Here 
we review the current understanding of the transition to jamming and the nature of the jammed 
state for disordered packings of particles that act through repulsive contact interactions and are at 
zero temperature and zero shear stress. We first discuss the breakdown of affine assumptions that 
underlies the rich mechanics near jamming. We then extensively discuss jamming of frictionless 
soft spheres. At the jamming point, these systems are marginally stable (isostatic) in the sense 
of constraint counting, and many geometric and mechanical properties scale with distance to this 
jamming point. Finally we discuss current explorations of jamming of frictional and non-spherical 
(ellipsoidal) particles. Both friction and asphericity tune the contact number at jamming away from 
the isostatic limit, but in opposite directions. This allows one to disentangle distance to jamming 
and distance to isostaticity. The picture that emerges is that most quantities are governed by the 
contact number and scale with distance to isostaticity, while the contact number itself scales with 
distance to jamming. 

PACS numbers: 61.43.-j, 64.70.D-, 83.80.Fg, 83.80.Hj, 83.80.Iz 



I. INTRODUCTION 

Jamming governs the transition to rigidity of disor- 
dered matter. Foams, emulsions, colloidal suspensions, 
pastes, granular media and glasses can jam in rigid, dis- 
ordered states in which they respond essentially elasti- 
cally to small applied shear stresses (Fig. QJi-d). How- 
ever, they can also easily be made to yield (unjam) and 
flow by tuning various control parameters. 

The transition from the freely flowing to the jammed 
state, the jamming transition, can be induced by varying 
thermodynamic variables, such as temperature or den- 
sity, but also mechanical variables such as the stress ap- 
plied to the sample: colloidal suspensions become col- 
loidal glasses as the density is increased near random 
close packing, flowing foams become static as the shear 
stress is decreased below the yield stress, and supercooled 
liquids form glasses as the temperature is lowered be- 
low the glass transition temperature. In 1998 Liu and 
Nagel presented their provocative jamming phase dia- 
gram (Fig. [TJ)), and proposed to probe the connections 
between various transitions to rigidity 

This review provides an overview of the current (par- 
tial) answers to the following two questions: What is the 
nature of the jammed state? What is the nature of the 
jamming transition? We focus on jammed model systems 
at zero temperature and zero shear — models for non- 
brownian emulsions, foams and granular media rather 
than colloidal and molecular glasses — and review the 
geometrical and mechanical properties of these systems 
as a function of the distance to jamming. 

In view of the very rapid developments in the field, 
the paper focuses on the basic jamming scenarios, which 
arises in (weakly) compressed systems of soft particles 



interacting through repulsive contact forces at zero tem- 
perature and zero shear. The picture that has emerged 
for the jamming transition in these systems is sufficiently 
complete to warrant an overview article and, in addition, 
provides a starting point for work on a wider range of 
phenomena, such as occurring in attractive systems 0], 
systems below jamming 14 , the flow of disordered media 
near jammin g [a , la, III , |8|, lilK jamming o f sy stems at finite 
temperature jlOl . Ill| and experiments [HI, El, EH ■ 

In this review the focus is on jamming of frictionless 
spheres, frictional spheres, and frictionless ellipsoids - 
soft (deformable) particles which interact through repul- 
sive contact forces. The distance to jamming of all these 
systems is set by the amount of deformation of the par- 
ticles, which can be controlled by the applied pressure 
or enforced packing fraction. These systems lose rigid- 
ity when the deformations vanish, or equivalcntly, when 
the confining pressure reaches zero. As we will see, these 
seemingly simple systems exhibit rich and beautiful be- 
havior, where geometry and mechanical response are in- 
tricately linked. 

The contact number, z, defined as the average num- 
ber of contacts per particle, plays a crucial role for these 
systems. There is a minimal value of z below which 
the system loses rigidity: when the contact number is 
too small, there are collective particle motions, so-called 
floppy modes, that (in lowest order) do not cost elastic 
energy. By a constraint counting argument one can es- 
tablish a precise value for the minimum value of z where 
the system does not generically allows floppy deforma- 
tions — this is the isostatic contact number Z{ so . As we 
will see, a host of mechanical and geometrical properties 
of jammed systems scale with distance to the isostatic 
point. 



2 




Temperature (e) 




1/Density 



FIG. 1: (a-d) Examples of everyday disordered media in a jammed state, (a) Granular media, consisting of solid grains in gas 
or vacuum, (b) Toothpaste, a dense packing of (colloidal) particles in fluid, (c) Mayonnaise, an emulsion consisting of a dense 
packing of (oil) droplets in an immiscible fluid, (d) Shaving foam, a dense packing of gas bubbles in fluid, (e) Jamming diagram 
proposed by Liu, Nagel and co-workers [J, Q] • The diagram illustrates that many disordered materials are in a jammed state 
for low temperature, low load and large density, but can yield and become unjammed when these parameters are varied. In 
this review we will focus on the zero temperature, zero load axis. For frictionless soft spheres, there is a well defined jamming 
transition indicated by point "J" on the inverse density axis, which exhibits similarities to an (unusual) critical phase transition. 



The crucial, and at first glance very puzzling point, is 
that while frictionless spheres reach isostaticity at the 
jamming point, frictional spheres are generally hyper- 
static {z > Zi so ) at jamming, while frictionless ellipsoids 
are hypostatic (z < Zi ao ) at jamming. As we will see, the 
relations between contact numbers, floppy modes, rigid- 
ity and jamming are subtle. 

Truly new and surprising physics emerges near jam- 
ming in systems as seemingly simple as disordered pack- 
ings of frictionless, deformable particles @. We first dis- 
cuss the breakdown of affinc assumptions that underlies 
the rich physics of jamming in section [TTJ We give an 
overview of the main characteristics of the jamming tran- 
sition for soft frictionless spheres in section IIIII Both 
friction and asphericity lead to new physics, as here the 
jamming transition and isostaticity decouple. Jamming 
of frictional soft spheres is discussed in section HVl and 
jamming of frictionless soft ellipsoids in section [V] Fi- 
nally, in section lVTl we sketch a number of open problems. 

II. MOTIVATION: MECHANICS OF 
DISORDERED MATTER 

The crucial question one faces when attempting to de- 
scribe the mechanics of materials such as foams, emul- 
sions or granular media, is how to deal with disorder. The 
simplest approach is to ignore disorder altogether, and at- 
tempt to gain insight based on models for ordered, "crys- 
talline" packings. A related approach, effective medium 
theory, does not strictly require ordered packings, but 



assumes that local deformations and forces scale sim- 
ilarly as global deformations and stresses. As we will 
see in section Til Al major discrepancies arise when these 
approaches are confronted with (numerical) experiments 
on disordered systems. This is because the response of 
disordered packings becomes increasingly non-afime near 
jamming (section III Bp . 



A. Failure of Affine Approaches 

1. Foams and Emulsions 

Some of the earliest studies that consider the ques- 
tion of rigidity of packings of particles concern the loss 
of rigidity in foams and emulsions with increasing wet- 
ness. Foams are dispersions of gas bubbles in liquid, 
stabilized by surfactant, and the gas fraction <f> plays 
a crucial role for the structure and rigidity of a foam. 
The interactions between bubbles are repulsive and vis- 
cous, and static foams are similar to the frictionless soft 
spheres discussed in section IIIII In real foams, gravity 
(which causes drainage) and gas diffusion (which causes 
coarsening) play a role, but we will ignore these. 

The unjamming scenario for foams is as follows. When 
the gas fraction approaches one, the foam is called dry. 
Application of deformations causes the liquid films to be 
stretched, and the increase in surface area then provides 
a restoring force: dry foams are jammed. When the gas 
fraction is lowered and the foam becomes wetter, the gas 
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bubbles become increasingly spherical, and the foam loses 
rigidity for some critical gas fraction <p c where the bubbles 
lose contact (Fig. [2J. The unjamming transition is thus 
governed by the gas fraction, which typically is seen as a 
material parameter. For emulsions, consisting of droplets 
of one fluid dispersed in a second fluid and stabilized by 
a surfactant, the same scenario arises. 

Analytical calculations are feasible for ordered pack- 
ings, because one only needs to consider a single particle 
and its neighbors to capture the packing geometry and 
mechanical response of the foam — due to the periodic 
nature of the packing, the response of the material is 
affine. The afhne assumption basically states that lo- 
cally, particles follow the globally applied deformation 
field — as if the particles are pinned to an elastically 
deforming sheet. More precisely, the strict definition of 
affine transformations states that three collincar particles 
remain collincar and that the ratio of their distances is 
preserved, and affine transformations are, apart from ro- 
tations and translations, composed of uniform shear and 
compression or dilatation. 

Packings of monodisperse bubbles in a two- 
dimensional hexagonal lattice ("liquid honeycomb" 
fl(| ) deform affincly. The bubbles lose contact at the 
critical density <j) c equal to w 0.9069, and ordered 

foam packings are jammed for larger densities [l6|, [l7| . 
When for such a model foam (f> is lowered towards 4> Cl 
the yield stress and shear modulus remain finite, and 
jump to zero precisely at d>r (T5, [l7T| . The contact number 
(average number of contacting neighbors per bubble) 
remains constant at 6 in the jammed regime. Similar 
results can be obtained for three-dimensional ordered 
foams, where 4> c is given by the packing density of the 
HCP lattice w 0.7405. 

3V2 

Early measurements for polydisperse emulsions by 
Princen and Kiss in 1985 [l8[ found a shear modulus 
which varied substantially with 4>. Even though no data 
was presented for <f> less than 0.75 and the fit only in- 
cluded points for which <f> > 0.8, the shear modulus was 
fitted as G ~ </> 1 / 3 (0 — cf> c ), where <p c « 0.71, and thus 
appeared to vanish at a critical density below the value 
predicted for ordered lattices [HI]. 

The fact that the critical packing density for ordered 
systems is higher than that for disordered systems may 
not be a surprise, given that at the jamming threshold, 
the particles are undeformed spheres, and it is well known 
that ordered sphere packings are denser than irregular 
ones However, the differences between the varia- 

tion of the moduli and yield strength with distance to 
the rigidity threshold predicted for ordered packings and 
measured for disordered emulsions strongly indicates that 
one has to go beyond models of ordered packings. 

2. Effective Medium Theory for Granular Media 

For granular media an important question has been 
to predict the bulk elasticity, and Makse and co-workers 
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FIG. 2: Simulated foam for increasing wetness, approaching 
unjamming for cf> [ 0.84 (adapted from pl| — Copyright by 
the American Physical Society). 

have carried out extensive studies of the variation of the 
elastic moduli and sound propagation speed with pres- 
sure in granular media from the perspective of effective 
medium theory [U EH, H2] • 

Effective medium theory (EMT) basically assumes 
that: (i) Macroscopic, averaged quantities can be ob- 
tained by a simple coarse graining procedure over the 
individual contacts, (ii) The effect of global forcing, e.g., 
imposing a deformation, trivially translates to changes in 
the local contacts. This second assumption is the "affine 
assumption" , and this will be the crucial assumption that 
breaks down near jamming. 

Makse et al. studied the breakdown of effective medium 
theory in the context of granular media. Assuming a 
Hertzian interaction between spherical grains (23[, the 
contact force / scales with the overlap 6 between particles 
as / ~ (5 3 / 2 . As a result, the stiffness of these contacts 
then scales as dsf ~ S 1 ^ 2 . Since, in good approximation, 
the pressure P ~ /, one obtains that the stiffness of the 
individual contacts scales as P 1 / 3 . EMT then predicts 
that the elastic bulk modulus K and shear modulus G 
scale as the stiffness of the contacts: K ~ G ~ P 1 / 3 , and 
that the sound velocities scales as P 1/6 [13, El E2, E3 . 
In particular, the ratio G/K should be independent of 
pressure. 

From a range of simulations Makse et al. concluded 
that the affine assumption works well for the compres- 
sion modulus provided that the change in contact num- 
ber with P is taken into account, but fails for the shear 
modulus — and suggested that this is due to the non- 
affine nature of the deformations [13, EH, E2|- We will 
discuss this issue at length in section IIIII 

B. Beyond AfHne Approaches 

In a seminal paper in 1990, Bolton and Weaire asked 
how a disordered foam loses rigidity when its gas fraction 
is decreased [ToT ]. They probed this question by simu- 
lations of a two-dimensional polydisperse foam, consist- 
ing of approximately hundred bubbles, as a function of 
<j> (Fig. [2]). Their model captures the essential surface 
tension driven structure of foams and predates the now 
widely used "surface evolver" code for foams [26j . 

The following crucial observations are made: (i) The 
critical density is around 0.84, which is identified as the 
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FIG. 3: Square root scaling of contact number z with 6 — <j) c 
observed in the Durian bubble model (adapted from [25( - 
Copyright by the American Physical Society). 



random close packing density in two dimensions — here 
the yield stress appears to vanish smoothly, (ii) The 
contact number z smoothly decreases with <f>. At (/) = 1 
the contact number equals six. This can be understood 
by combining Eulcr's theorem which relates the number 
of vertices, faces and edges in tilings with Plateau's rule 
that for a two-dimensional dry foam in equilibrium, three 
films (faces) meet in one point (vertex). When <j> — > <j) c , 
the contact number appears to reach the marginal value, 
four, (in) The shear modulus decreases with <f> and ap- 
pears to smoothly go to zero at <j> = <f> c (unfortunately 
the authors do not comment on the bulk modulus). 

In related work on the so-called bubble model de- 
veloped for wet foams in 1995, Durian reached similar 
conclusions for two-dimensional model foams, and more- 
over found that the contact number indeed approaches 
4(= 2d) near jamming, and observed the non-trivial 
square root scaling of z — 4 with excess density for the 
first time (FigEJ). All these findings are consistent with 
what is found in closely related models of frictionless soft 
spheres near jamming, as discussed in sections Mil 

Experimentally, measurements of the shear modulus 
and osmotic pressure of compressed three-dimensional 
monodisperse but disordered emulsions found similar be- 
havior for the loss of rigidity [27|, [H, [2{|. The shear 
modulus, (when scaled appropriately with the Laplace 
pressure, which sets the local "stiffness" of the droplets) 
grows continuously with and vanishes at <p c ps 0.635, 
corresponding to random close packing in three dimen- 
sions. The osmotic pressure exhibits very similar scaling, 
implying that the bulk modulus (being proportional to 
the derivative of the pressure with respect to <fi) scales 
differently from the shear modulus — the difference be- 
tween shear and bulk modulus is another hallmark of 
jamming of frictionless spheres. 




FIG. 4: Deformation fields of packings of 1000 frictionless par- 
ticles under compression (a,c) and shear (b,d) as indicated by 
the red arrows. The packings in the top row (a,b) are strongly 
jammed (contact number z = 5.87), while the packings in the 
bottom row (c,d) are close to the jamming point — their con- 
tact number is 4.09, while the jamming transition occurs for 
z — 4 in this case. Clearly, the deformation field becomes in- 
creasingly non-affine when the jamming point is approached 
(adapted from [30l.l3lT| — Copyright by the American Physical 
Society). 

There is thus a wealth of simulational and experimen- 
tal evidence that invalidates simple predictions for the 
rigidity of disordered media based on our intuition for 
ordered packings. The crucial ingredient that is missing 
is the non-affine nature of the deformations of disordered 
packings (Fig. . There is no simple way to estimate the 
particles motion and deformations in disordered systems, 
and one needs to resort to (numerical) experiments. Jam- 
ming can be seen as the avenue that connects the results 
of such experiments. Jamming aims at capturing the me- 
chanical and geometric properties of disordered systems, 
building on two insights: first, that the non-affine char- 
acter becomes large near the jamming transition, and 
second, that disorder and non-affinity are not weak per- 
turbations away from the ordered, affine case, but may 
lead to completely new physics f^lHEllHElHllaEEHl- 

III. JAMMING OF SOFT FRICTIONLESS 
SPHERES 

Over the last decade, tremendous progress has been 
made in our understanding of what might be considered 
the "Ising model" for jamming: static packings of soft, 
frictionless spheres that act through purely repulsive con- 
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FIG. 5: States of soft frictionless spheres as function of pack- 
ing density (j>, below, at, and above the critical density <j> c . 
Left: Unjammed system at a density below the critical den- 
sity — pressure is zero and there are no contacts. Middle: 
Marginally rigid system consisting of undeformed frictionless 
spheres just touching. The system is at the jamming tran- 
sition (point J), has vanishing pressure, critical density and 
2d contacts per particle, where d is the dimension. Right: 
Jammed system for finite pressure and density above (j> c - 



tact forces. In this model, temperature, gravity and shear 
are set to zero. The beauty of such systems is that they 
allow for a precise study of a jamming transition. As we 
will see in sections |IV] and [V] caution should be applied 
when applying the results for soft frictionless spheres to 
frictional and/or non-spherical particles. 

From a theoretical point of view, packings of soft fric- 
tionless spheres are ideal for three reasons. First, they 
exhibit a well defined jamming point: For positive P the 
system is jammed, as it exhibits a finite shear modulus 
and a finite yield stress 0] , while at zero pressure the sys- 
tems loses rigidity. Hence, the (un)jamming transition 
occurs when the pressure P approaches zero, or, geomet- 
rically, when the deformations of the particles vanish. 
The zero pressure, zero shear, zero temperature point in 
the jamming phase diagram is referred to as "point J" 
(Fig. [TJ? and [5]). In this review, point J will only refer 
to soft frictionless spheres and not to jamming transi- 
tions of other types of particles. Second, at point J the 
contact number approaches the so-called isostatic value, 
and the system is marginally stable. The system's me- 
chanical and geometrical properties are rich and peculiar 
here. For large systems the critical packing density, </> c , 
approaches values usually associated with random close 
packing. Third, the mechanical and geometrical proper- 
ties of jammed systems at finite pressure, or cquivalcntly, 
4> — 4> c > 0, exhibit non-trivial power law scalings as a 
function Acjj := <f> — <f> c or, similarly, as function of the 
pressure, P. 

In this section we address the special nature of point J 
and discuss the scaling of the mechanical and geometrical 
properties for jammed systems near point J. We start in 
section IIII Al by a brief discussion of a few common con- 
tact laws and various numerical protocols used to gener- 
ate jammed packings. We then present evidence that the 
jamming transition of frictionless spheres is sharp and 
discuss the relevant control parameters in section IIII Bl 



In section IIII CI we discuss the special geometrical fea- 
tures of systems at point J, as probed by the contact 
number and pair correlation function. Away from point 
J the contact number exhibits non-trivial scaling, which 
appears to be closely related to the pair correlation func- 
tion at point J, as discussed in section UlIDI Many fea- 
tures of systems near point J can be probed in linear re- 
sponse, and these are discussed at length in section Ull El 
— these include the density of states (|III E 1[) , diverging 
length and time scales (jlll E 2p . elastic moduli (jlll E 3p 
and non-affine displacements (|III E 4[) . We close this sec- 
tion by a comparison of effective medium theory, rigidity 
percolation and jamming, highlighting the unique nature 
of jamming near point J piIE 5|) . 

A. Definition of the Model 

At the (un)jamming transition soft particles are un- 
deformed, and the distance to jamming depends on the 
amount of deformation. Rigid particles are therefore al- 
ways at the jamming transition, and soft particles are 
necessary to vary the distance to point J. Deformable fric- 
tionless spheres interact through purely repulsive body 
centered forces, which can be written as a function of 
the amount of virtual overlap between two particles in 
contact. Denoting the radii of particles in contact as Ri 
and Rj and the center-to-center distance as ry , it is con- 
venient to define a dimcnsionlcss overlap parameter Sij 
as 



so that particles are in contact only if Sij > 0. We limit 
ourselves here to interaction potentials of the form: 

Vij = etj 5% S l3 > , (2) 

Vij = < . (3) 

By varying the exponent, a, one can probe the nature and 
robustness of the various scaling laws discussed below. 
For harmonic interactions, a = 2 and sets the spring 
constant of the contacts. Hertzian interactions between 
three-dimensional spheres, where contacts are stiffer as 
they are more compressed, correspond to a = 5/2 [9lj . 
O'Hcrn et al have also studied the "Herman" interaction 
(a = 3/2), which corresponds to contacts that become 
progressively weaker when compressed 0] ■ 

Once the contact laws are given, one can generate pack- 
ings by various different protocols, of which MD (Molec- 
ular Dynamics) [2fj|. l2ll. [23 . l2~i a nd conjugate gradient Q 
are the most commonly used \92\ . In MD simulations one 
typically starts simulations with a loose gas of particles, 
which are incrementally compressed, either by shrinking 
their container or by inflating their radii. Supplement- 
ing the contact laws with dissipation (inelastic collisions, 
viscous drag with a virtual background fluid, etc) the 
systems "cools" and eventually one obtains a stationary 
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jammed state. While straightforward, one might worry 
that statistical properties of packings obtained by such 
procedure depend on aspects of the procedure itself — 
for frictional packings, this is certainly the case [37| . 

For frictionless particles, the interactions are conser- 
vative, and one can exploit the fact that stable packings 
correspond to minima of the elastic energy. Packings 
can then be created by starting from a completely ran- 
dom configuration and then bringing the system to the 
nearest minimum of the potential energy. When the en- 
ergy at this minimum is finite, the packing is at finite 
pressure, and this procedure is purported to sample the 
phase space of allowed packings flatly 0, HI]. An ef- 
fective algorithm to find such minima is known as the 
"conjugate gradient technique" (39|. For frictionless sys- 
tems, we are not aware of significant differences between 
packings obtained by MD and by this method [93| . 

Finally it should be noted that to avoid crystallization, 
two-dimensional packings are usually made polydispcrse, 
and a popular choice are bidisperse packings where parti- 
cles of radii 1 and 1.4 are mixed in equal amounts @,H3]- 
In three dimensions, this is not necessary as monodis- 
perse spheres then do not appear to order or crystal- 
lize for typically employed numerical packing generation 
techniques. 



B. Evidence for Sharp Transition 

The seminal work of O'Hern et al. 0, H(| has laid the 
groundwork for much of what we understand about jam- 
ming of frictionless soft spheres. These authors begin 
by carefully establishing that frictionless soft spheres ex- 
hibit a sharp jamming transition. First, it was found 
that when a jammed packing is decompressed, the pres- 
sure, the bulk modulus and the shear modulus vanish 
at the same critical density <fi c . For finite systems, the 
value of 4>c varies from system to system. For systems of 
1000 particles the width of the distribution of <fi c , W, still 
corresponds to 0.4%, and must therefore not be ignored. 
Second, it was shown that the width, W, vanishes with 
the number of particles N as W ~ N -1 / 2 — independent 
of dimension, interaction potential or polydispersity. In 
addition, the location of the peak of the distribution of 
</> c , </>o, also scales with N: fo-c/)* = (0.12±0.03)7V" 1 / i/d . 
Here d is the dimensionality, v = 0.71 ± 0.08 and <f>* 
approaches 0.639 ± 0.001 for three-dimensional monodis- 
perse systems. 

These various scaling laws suggest that for friction- 
less spheres the jamming transition is sharp in the limit 
of large systems. This jamming point is referred to as 
point J (see Fig. and [5]). At the jamming point, the 
packings consist of perfectly spherical (i.e., undeformed) 
spheres which just touch (Fig. [5]). The packing fraction 
for large systems, <fi* , reaches values which have been 
associated with random close packing (RCP) [3, [l5| - 
(~ 0.84 in two dimensions, ~ 0.64 in three dimensions). 
It should be noted that the RCP concept itself is contro- 



versial fill ]. 

Control Parameters — As we will see, the properties 
of packings of soft slippery balls are controlled by their 
distance to point J. What is a good control parameter 
for jamming at point J? The spread in critical density 
for finite systems indicates that one should not use the 
density, but only the excess density Acj> := <j> — <j> c as 
control parameter. In other words, fixing the volume is 
not the same as fixing the pressure for finite systems. 

The disadvantage of using the excess density is that 
it requires deflating packings to first obtain <p c This 
extra step is not necessary when P is used as control 
parameter, since the jamming point corresponds to P = 
no matter what the system size or <j) c is of a given 
system. While we believe it is much simpler to deal with 
fixed pressure than with fixed volume, a disadvantage of 
P is that its relation to A0 is interaction dependent: the 
use of the excess density stresses the geometric nature of 
the jamming transition at point J. 

We suggest that the average overlap (8) is the simplest 
control parameter — even though its use is not common. 
First, (5) is geometric and interaction independent and 
reaches zero at jamming, also for finite systems. More- 
over, for finite systems (S) still controls the pressure and 
will be very close to A(j). Of course, in infinite systems, 
control parameters like the pressure P, the average parti- 
cle overlap (5) and the density <f> are directly linked — for 
interactions of the form Eq. ©, P ~ S"' 1 ~ (A^)"" 1 . 
Below, we will use a combination of all these control pa- 
rameters, reflecting the different choices currently made 
in the field. 



C. Geometry at Point J 

At point J, the system's packing geometry is hi ghly 
non-trivial. First, systems at point J are isostatic [43|: 
the average number of contacts per particle is sharply 
defined and equals the minimum required for stability 
0, 0, Second, near jammin g g (r) diverges when 

r | 1 (for particles of radius 1) [42LI46I. Ul\ . 

Isostaticity — The fact that the contact number at 
point J attains a sharply defined value has been argued to 
follow directly from counting the degrees of freedom and 
constraints |4f| . We discuss such counting arguments 
in detail in Appendix \X\ but give here the gist of the 
argument for frictionless spheres. 

Suppose we have a packing of N soft spheres in d 
dimensions, and that the contact number, the average 
number of contacts at a particle, equals z — the to- 
tal number of contacts equals zN/2, since each contact 
is shared by two particles. First, the resulting packing 
should not have any floppy modes, deformation modes 
that cost zero energy in lowest order. As we discuss in 
Appendix[Al this is equivalent to requiring that the Nz/2 
contact forces balance on all grains, which yields dN con- 
straints on Nz/2 force degrees of freedom: hence z > 2d. 
The minimum value of z required is referred to as the 
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FIG. 6: The pair correlation function g(r > 1) of a three- 
dimensional system of monodisperse spheres of radius 1, il- 
lustrates the abundance of near contacts close to jamming 
(A(f> = 10 -8 here). From [Z^] — Copyright by the American 
Physical Society. 




FIG. 7: (a) Excess contact number z — z c as function of ex- 
cess density cj> — 4> c . Upper curves: represent monodisperse 
and bidisperse packings of 512 soft spheres in three dimen- 
sions with various interaction potentials, while lower curves 
correspond to bidisperse packings of 1024 soft discs in two 
dimensions. The straight lines have slope 0.5. From || - 
Copyright by the American Physical Society, (b) Schematic 
contact number as function of density, illustrating the mixed 
nature of the jamming transition for frictionless soft spheres. 



isostatic value z; so : for frictionless spheres, Zj SO = 2d. 

Second, at point J, since the particles are undcformcd: 
the distance between contacting particles has to be pre- 
cisely equal to the sum of their radii. This yields Nz/2 
constraints for the dN positional degrees of freedom: 
therefore, one only expects generic solutions at jamming 
when z < 2d. 

Combining these two inequalities then yields that the 
contact number z c at the jamming point for soft fric- 
tionless disks generically will attain the isostatic value: 
z c = z iso = 2d 0, H, El. As we will see below, such 
counting arguments should be regarded with caution, 
since they do not provide a correct estimate for the con- 
tact number at jamming of frictionless ellipsoidal parti- 
cles [II, El, [Hoi. 

Numerically, it is far from trivial to obtain convincing 
evidence for the approach of the contact number to the 
isostatic value. Apart from corrections due to finite sys- 
tem sizes and finite pressures, a subtle issue is how to 
deal with rattlers, particles that do not have any con- 
tacts with substantial forces, but still arise in a typical 
simulation. These particles have low coordination num- 
ber and their overlap with other particles is set by the 
numerical precision — these particles do not contribute 
to rigidity. For low pressures, they can easily make up 
5% of the particles. An accurate estimate of the contact 
number than requires one to ignore these particles and 
the corresponding "numerical" contacts pi. FtoI] . 

Pair Correlation Function In simulations of 

monodisperse spheres in three dimensions, it was found 
that near jamming g(r) diverges when r [ \ (for particles 
of radius 1 ) : 



This expresses that at jamming a singularly large number 



of particles are on the verge of making contact (Fig. [5]) 
[43 . [46[ . This divergence has also been seen in pure hard 
sphere packings [47fl , In addition to this divergence, g(r) 
exhibits a delta peak at r = 1 corresponding to the dN/2 
contacting pairs of particles. 

In simulations of two-dimensional bidisperse systems, a 
similar divergence can be observed, provided one studies 
.9(0' where the rcscalcd intcrparticlc distance £ is defined 
as r I (Ri + Rj), and where Ri and Rj are the radii of the 
undeformed particles in contact [5l| . 



D. Relating Contact Numbers and Packing 
Densities away from J 

Below jamming, there are no load bearing contacts and 
the contact number is zero, while at point J, the contact 
number attains the value 2d. How does the contact num- 
ber grow for systems at finite pressure? Assuming that 
(i) compression of packings near point J leads to essen- 
tially affine deformations, and that (ii) g(r) is regular 
for r > 1, z would be expected to grow linearly with <p: 
compression by 1% would then bring particles that are 
separated by less than 1% of their diameter in contact, 
etc. But we have seen above that g(r) is not regular, and 
we will show below that deformations are very far from 
affine near jamming — so how does z grow with </>? 

Many authors have found that the contact number 
grows with the square root of the excess density A</> := 
(j)-(j)c HUES!!! ( see Fi g- 01 ■ O'Hern et al. have stud- 
ied this scaling in detail, and find that the excess contact 
number Az := z — z c scales as Az ~ (/\0)O.5O±o.o3^ wnere 
z c , the critical contact number, is within error bars equal 
to the isostatic value 2d Q. Note that this result is in- 
dependent of dimension, interaction potential or polydis- 
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persity (see Fig.[7ji). Hence, the crucial scaling law is 

Az = z 0V fA4>, (5) 

where the precise value of the prcfactor zq depends on 
dimension, and possibly weakly on the degree of polydis- 
persity, and is similar to 3.5 ± 0.3 in two dimensions and 
7.9 ± 0.5 in three dimensions Q. 

The variation of the contact number near J can there- 
fore be perceived to be of mixed first/second order char- 
acter: below jamming z = 0, at J the contact number z 
jumps discontinuously from zero to 2d, and for jammed 
systems the contact number exhibits non-trivial power 
law scaling as a function of increasing density (Fig. [3] 
and El). 

We will see below that many other scaling relations 
(for elastic moduli, for the density of state and for char- 
acteristic scales) are intimately related to the scaling of 
z, and the contact number scaling can be seen as the cen- 
tral non-trivial scaling in this system. (In frictional and 
non spherical packings, similar scalings for z are found.) 

A subtle point is that the clean scaling laws for Az vs 
A<f> are only obtained if one excludes the rattlers when 
counting contacts, but includes them for the packing frac- 
tion 0. Moreover, for individual packings the scatter in 
contact numbers at given pressure is quite substantial — 
see for example Fig. 9 from [52j — and smooth curves 
such as shown in Fig. [7^ can only obtained by averag- 
ing over many packings. Finally, the density <p is usually 
defined by dividing the volume of the undeformed parti- 
cles by the box size, and packing fractions larger than 1 
are perfectly reasonable. Hence, in comparison to pack- 
ing fractions defined by dividing the volume of the de- 
formed particles by the box size, <p is larger because the 
overlap is essentially counted double. Even though none 
of these subtleties should play a role for the asymptotic 
scaling close to jamming in large enough systems, they 
are crucial when comparing to experiments and also for 
numerical simulations. 



1. Connections between contact number scaling, g(r) and 
marginal stability 

The scaling of Az can be related to the divergence of 
the radial distribution function as follows [56| . Imagine 
compressing the packing, starting from the critical state 
at point J, and increasing the typical particle overlap 
from zero to <5. If one assumes that this compression 
is essentially affine, then it is reasonable to expect that 
such compression closes all gaps between particles that 
are smaller than 5. Hence 

Wyart approaches the square root scaling of Az from a 
different angle, by first showing that the scaling Az ~ V5 
is consistent with the system staying marginally stable 




FIG. 8: Definition of relative displacement Uij, uu and u±. 



at all densities, and then arguing that the divergence 
in g(r) is a necessary consequence of that [54[. Both 
his arguments require assumptions which are not self- 
evident, though j52[. 



E. Linear Response and Dynamical Matrix 

A major consequence of isostaticity at point J is that 
packings of soft frictionless spheres exhibit increasingly 
anomalous behavior as the jamming transition is ap- 
proached. That anomalies occur near jamming is ulti- 
mately a consequence of the fact that the mechanical re- 
sponse of an isostatic system cannot be described by elas- 
ticity — isostatic systems are essentially different from 
ordinary elastic systems [Ifl HH . 

In principle these anomalies can be studied at the jam- 
ming point, however, much insight can be gained by ex- 
ploring the mechanical properties as a function of dis- 
tance to the isostatic point. Below we review a number of 
such non-trivial behaviors and scaling laws that arise near 
point J. We will focus on the response to weak quasistatic 
perturbations, and on the vibrational eigenfrequencies 
and eigenmodes of weakly jammed systems. Both are 
governed by the dynamical matrix of the jammed pack- 
ing under consideration. 

For linear deformations, the changes in elastic energy 
can be expressed in the relative displacement u.y of neigh- 
boring particles i and j. It is convenient to decompose 
Uij in components parallel (mm) and perpendicular 
to Tij, where connects the centers of particles i and 
j (Fig. [8]). In these terms the change in energy takes a 
simple form HH, S3, , 

where and fey denote the contact forces and stiff- 
nesses. For power law interactions of the form given in 
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FIG. 9: Density of vibrational states D(u) for 1024 spheres 
interacting with repulsive harmonic potentials. Distance to 
jamming A(f> equals 0.1 (black), 10~ 2 (blue), 10~ 3 (green), 
10 -4 (red) and 10 -8 (black). The inset shows that the char- 
acteristic frequency cj*, defined as where D(ui) is half of the 
plateau value, scales linearly with Az. The line has slope 
1. Adapted from [Hi], [5jJ — Copyright by the Institute of 
Physics 

Eq. j2]), wc can rewrite this as [3(ij |: 

A£ = ^^(^-^) ' (8) 

The dynamical matrix Aiij^p is obtained by rewriting 
cq. ([7]) in terms of the independent variables, Ui.„, as 

Here M. is a dN x dN matrix with N the number of 
particles, indices n,m label the coordinate axes, and the 
summation convention is used. 

The dynamical matrix contains all information on the 
clastic properties of the system. By diagonalizing the 
dynamical matrix one can probe the vibrational proper- 
ties of systems near jamming 0, [33|, H3, [Ei| (see section 
IIIIE 1[) . The dynamical matrix also governs the elastic 
response of the system to external forces / oxt (see sec- 
tions mminiEi) \m,\Ej§-- 

•Mij,nm ^j,m = fin ' (~^) 

1. Density of States 

Studies of the vibrational modes, and the associated 
density of (vibrational) states (DOS) have played a key 
role in identifying anomalous behavior near point J. Low 
frequency vibrations in ordinary crystalline or amorphous 
matter are long-wavelength plane waves. Counting the 
number of these, one finds that the density of vibrational 
states D{lo) is expected to scale as D{ui) ~ for low 
frequencies — this is called Debye behavior. Jammed 
packings of frictionless spheres do show Debye-like be- 
havior far away from jamming, but as the point J is ap- 
proached, both the structure of the modes and the den- 
sity of states exhibit surprising features @, H3, HE HI] . 



The most striking features of the density of states are 
illustrated in Fig. [HI First, far above jamming, the DOS 
for small frequencies is regular (black curve). Second, ap- 
proaching point J, the density of vibrational states DOS 
at low frequencies is strongly enhanced. (In analogy to 
what is observed in glasses, this is sometimes referred 
to as the boson peak, since the ratio of the observed 
DOS and the Debye prediction exhibits a peak at low to). 
More precisely, the DOS becomes essentially constant up 
to some low-frequency crossover scale at uo = w* , be- 
low which the continuum scaling ~ w^ 1 is recovered. 
Third, the characteristic frequency lo* vanishes at point 
J as to* ~ Az. 

The density of states thus convincingly shows that, 
close to the isostatic point / jamming point, the material 
is anomalous in that it exhibits an excess of low frequency 
modes, and that at point J, the material does not ap- 
pear to exhibit any ordinary Debye/continuum behavior 
as here the DOS becomes flat. Jamming of frictionless 
spheres thus describes truly new physics. 

Normal Modes — The nature of the vibrational modes 
changes strongly with frequency, and, to a lesser extend, 
with distance to point J. Various order parameters can 
be used to characterize these modes, such as the (inverse) 
participation ratio, level repulsion and localization length 
[58l . [59| . The participation ratio for a given mode is de- 
fined as V = {1/N) (T,i\ Ul \ 2 ) 2 /_Ei\ Ui \ 4 , where u % is the 
polarization vector of particle i [581 ] . It characterizes how 
evenly the particles participate in a certain vibrational 
mode — extended modes have V of order one, while lo- 
calized modes have smaller V 1 with hypothetical modes 
where only one particle participates reaching V = 1/N. 

Studies of such order parameters have not found very 
sharp changes in the nature of the modes either with 
distance to jamming or with eigenfrequency [58l. IBl^. |60| . 
It appears to be more appropriate to think in terms of 
typical modes and crossovers. Qualitatively, one can con- 
sider the DOS to consist of roughly three bands: a low 
frequency band where D{uj) ~ u> d ~ 1 , a middle frequency 
band where D(u) is approximately fiat, and a high fre- 
quency band where D(lu) decreases with u> [58j . 

Representative examples of modes in these three bands 
are shown in Fig. 1101 The modes in the low frequency 
band come in two flavors: plane wave like with V ~ 1, 
and quasi localized with small V [5^, [6(| ■ The modes in 
the large frequency band arc essentially localized with 
small V . The vast majority of the modes are in the 
mid frequency band (especially close to jamming), and 
are extended but not simple plain waves — typically the 
eigenvectors have a swirly appearance. 

The localization length £ of these modes has been es- 
timated to be large, so that many modes have £ com- 
parable or larger than the system size. Consistent with 
this, the modes in the low and mid frequency range are 
mostly extended, £ > L, and exhibit level repulsion (i.e., 
the level spacing statistics P(Auj) follows the so-called 
Wigncr surmise of random matrix theory), while the high 
frequency modes are localized (£ < L) and exhibit Pois- 
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FIG. 10: Representative eigenmodes for a two-dimensional system of 10 4 particles interacting with three-dimensional Hertzian 
interactions (a = 5/2, see Eq.[2]) at a pressure far away from jamming (z w 5.09). For all modes, the length of the vectors oc m 
are normalized such that <7i|tti| 2 is a constant, (a) Continuum-like low frequency mode at lu « 0.030, V « 0.79, and i u = 3, 
where i w counts the non trivial modes, ordered by frequency, (b) Quasi-localized low frequency mode at lu w 0.040, V ~ 0.06, 
and i u = 7. (c) Disordered, "swirly" mid frequency mode at lu m 0.39, V ~ 0.31, and i u = 1000. (d) Localized high frequency 
mode at lu « 4.00, P w 0.0013, and i u = 9970. 



sonian level statistics [59j . 

When point J is approached, the main change is that 
the low frequency, "Debye" range shrinks, and that both 
the number of plane waves and of quasi-localized reso- 
nances diminishes [58|, [59|, 16CJ • 

2. Characteristic Length and Time Scales 

The vanishing of the characteristic frequency lu* at 
point J suggests to search for a diverging length scale. 
Below we give an analytical estimate for this length scale 
and discuss indirect and direct observations of this length 
scale in simulations. 

Estimate of I* — As pointed out by Wyart et al. [54| . 
if we cut a circular blob of radius I from a rigid mate- 
rial, it should remain rigid. The rigidity (given by the 
shear modulus) of jammed materials is proportional to 
Az. The circular blob has of the order £ d Az excess con- 
tacts. By cutting it out, one breaks the contacts at the 
perimeter, of which there are of order z£ d ~ x . If the num- 
ber of broken contacts at the edge is larger than the num- 
ber of excess contacts in the bulk, the resulting blob is 
not rigid, but floppy: it can be deformed without energy 
cost (in lowest order). The smallest blob one can cut out 
without it being floppy is obtained when these numbers 
are equal, which implies that it has radius £* ~ zj Az. 
Close to the jamming transition, z is essentially constant, 
and so one obtains as scaling relation that [5J] 



Observation of I* in Vibration Modes — Using the 
speed of sound one can translate the crossover frequency 
lu* into a wavelength, which scales as At ~ 1/%/ Az for 
transverse (shear) waves and as Al ~ 1/ Az for longitudi- 
nal (compressional) waves — the difference in scaling is 
due to the difference in scaling of shear and bulk moduli 
(see section lTlI E 3l belowV By examining the spatial vari- 
ation of the eigenmodc corresponding to the frequency 
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FIG. 11: Divergence of a characteristic length scale near jam- 
ming as observed in the fluctuations of the changes of contact 
forces of a system of 10 4 Hertzian discs. Blue (red) bonds cor- 
respond to increased (decreased) force in response to pushing 
a single particle in the center of the packing to the right. In 
panel (a), the system is far from jamming and z = 5.55, while 
in panel (b), the system is close to jamming and z — 4.05 
(adapted from [3l|). 

lu* , At has been observed by Silbert et al. [56|. Notice, 
however, that the scaling of At is different from the scal- 
ing of I* — it is Al that coincides with the length scale 
I* derived above. 

Observation of I* in Point Response — The signature 
of the length scale £* can be observed directly in the point 
force response networks : Close to point J, i.e. for small 
Az, the scale up to which the response looks disordered 
becomes large (see Fig. [TT]) [3(1 EJ. By studying the 
radial decay of fluctuations in the response to an inflation 
of a single central particle (which is more symmetric than 
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particles i and j via the Irving-Kirkwood equation: 
1 
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FIG. 12: Bulk (K) and shear (G) modulus as function of 
distance to jamming for two-dimensional bidisperse systems, 
with interaction potential V ~ 5 a (see Eqs. [2}. The closed 
symbols denote moduli calculated by forcing the particles to 
move afnnely, and the open symbols correspond to the moduli 
calculated after the system has relaxed. Slopes as indicated 
(adapted from 0] — Copyright by the American Physical 
Society) . 



that of point forcing as shown in Fig. Ill)) as a function 
of distance to jamming, one obtains a crossover length I* 
which, as the theoretically derived length scale, varies as 
I* w 6/Az [H|. 

Characteristic Length and Validity of Elasticity — An 
important issue, which has in particular been studied ex- 
tensively in the context of granular media, is whether 
elasticity can describe a systems response to, for example, 
point forcing [HBI. f6ll| . Extensive observations of the linear 
response, connected to the direct observation of I*, sug- 
gest that there is a simple answer, and that the distance 
to the isostatic limit is crucial [lo, HH : Below a length 
scale I* the response is dominated by fluctuations, and 
the deformation field can be seen as a distorted floppy 
mode, while at larger length scales the systems response 
crosses over to elasticity. This is for a single realization 
— it can also be shown that, even close to jamming, the 
ensemble averaged response of a weakly jammed system 
is consistent with elasticity, provided the correct values 
of the elastic moduli are chosen — these moduli are con- 
sistent with the globally defined ones [3~D ]. 



3. Scaling of Shear and Bulk Moduli 

The scaling of the shear modulus, G, and bulk modu- 
lus, K, plays a central role in connecting the non-affinc, 
disordered nature of the response to the anomalous elas- 
tic properties of systems near jamming. To understand 
why disorder is so crucial for the global, mechanical re- 
sponse of collection of particles that act through short 
range interactions, consider the local motion of a packing 
of spherical, soft frictionless spheres under global forcing. 
The global stresses can be obtained from the relative po- 
sitions fij and contact forces /y of pairs of contacting 



(12) 



where <j a b is the stress tensor, a and (3 label coordinates, 
and V is the volume. 

Once we know the local motion of the particles in re- 
sponse to an externally applied deformation, we can cal- 
culate the contact forces from the force law and obtain 
thus the stress in response to deformation. Let us first 
estimate the scaling of the moduli from the affine predic- 
tion where one assumes that the typical particle overlap 
6 is proportional to A<f> and that all bonds contribute 
similarly to the increase in elastic energy when the pack- 
ing is deformed. For a deformation strain e we can esti- 
mate the corresponding increase in energy from Eq. @ 
as AE ~ Y^ke 2 . Therefore, under affine deformations, the 
corresponding elastic modulus is of order k — in other 
words, the elastic moduli simply follow from the typical 
stiffnesses of the contacts. 

Consider now deforming a disordered jammed pack- 
ing. All particles feel a local disordered environment, 
and deformations will not be affine (Fig. Q|. The point is 
that these non-affine motions become increasingly strong 
near the jamming transition, and qualitatively change 
the scaling behavior of, e.g., t he shear modulus of foams 
and granular media 0, Il5l. l20l l43l . [62j . 

A particularly enlightening manner to illustrate the 
role of non-affine deformations is to initially force the 
particle displacements to be affine, and then let them re- 
lax. In general, the system can lower its elastic energy 
by additional non-affine motions. Calculating the clas- 
tic energies of enforced affine deformations and of the 
subsequent relaxed packings of soft frictionless spheres, 
O'Hern and co-workers found that the non-affine relax- 
ation lowers both the shear and bulk modulus, but cru- 
cially changes the scaling of the shear modulus with dis- 
tance to jamming Q — see Fig. [TJ] 

In general, one finds that for power law interactions 
(Eq. [2]), the pressure scales as A</)" -1 and the contact 
stifness k and bulk modulus K scale as A<fr a ~ 2 pi. l30l. |62| . 
The surprise is that the shear modulus G gets progres- 
sively smaller than the bulk modulus near point J, and 
G scales differently from K with distance to jamming: 
G ~ A0"" 3 / 2 (See Fig.rjJ @, M, M, HI . The relations 
between the scaling of G, K and k can be rewritten as 



G - AzK ~ Az k 



(13) 



It is worth noting that many soft matter systems (pastes, 
emulsions) have shear moduli which arc much smaller 
than compressional moduli — from an application point 
of view, this is a crucial property. 

Putting all this together, we conclude that the affine 
assumption gives the correct prediction for the bulk mod- 
ulus (since k ~ 5 a ~ 2 ~ A</>" -2 ), but fails for the shear 
modulus. This failure is due to the strongly non-affinc 
nature of shear deformations: deviations from affine de- 
formations set the elastic constants 0, HO, HH, [62j]. 
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FIG. 13: (a) Illustration of definition of displacement angle a. (b-c) Probability distributions P(a) for compression (b) and 
shear (c), for Hertzian particles in two dimensions. The three pressures indicated correspond to z as 6.0, z as 4.5 and z as 4.1 
respectively (adapted from [3(| — Copyright by the American Physical Society). 



As we will see below, the correspondence between the 
bulk modulus and the affine prediction is fortuitous, since 
the response becomes singularly non-affine close to point 
J, for both compressive and shear deformations (section 
IIII E 51) . 



4- Non-Affine Character of Deformations 

Approaching the jamming transition, the spatial struc- 
ture of the mechanical response becomes less and less 
similar to continuum elasticity, but instead increas- 
ingly reflects the details of the underlying disordered 
packing and becomes increasingly non-affine [30| — see 
Fig. [3Jt. Here we will discuss this in the light of 
Eq. ©, which expresses the changes in energy as func- 
tion of the local deformations uii and u± : AE = 

To capture the degree of non-affinity of the response, 
Ellcnbrock and co-workers have introduced the displace- 
ment angle onj [94j . Here a^- denotes the angle between 
Uii and rii, or, 



tancvij = 



(14) 



The probability distribution P(a) can probe the degree 
of non-affinity by comparison with the expected P(a) for 
affine deformations. Affine compression corresponds to 
a uniform shrinking of the bond vectors, i.e. u±jj = 
while tin ij = —£r.ij < 0: the corresponding P(a) exhibits 
a delta peak at a — it. The effect of an affine shear on a 
bond vector depends on its orientation, and for isotropic 
random packings, P(a) is flat. 

Numerical determination of P(ot) show that systems 
far away from the jamming point exhibit a P(a) similar 
to the affine prediction, but that as point J is approached, 
P(a) becomes increasingly peaked around a = ir/2 



(Fig. (|T3|) b-c). This is reminiscent of the P(a) of floppy 
deformations, where the bond length does not change 
and P(a) exhibits a 5-peak at ir/2. Hence deformations 
near jamming become strongly non-affine, and, at least 
locally, resemble those of floppy modes. 

Non-affinity of Floppy Modes and Elastic Response — 
Wyart and co-workers have given variational arguments 
for deriving bounds on the energies and local deforma- 
tions of soft (low energy) modes starting from purely 
floppy (zero energy) modes [H, [63[ . They construct trial 
soft modes that are basically floppy modes, obtained by 
cutting bonds around a patch of size I* , and then mod- 
ulating these trial modes with a sine function of wave- 
length I* to make the displacements vanish at the loca- 
tions of the cut bonds [30l . 1 5 411 . In particular, for the 
local deformations, they find [63| 



"II 
u± 



1 



"II 
u± 



Az, 



(15) 



where symbols without indices ij refer to typical or av- 
erage values of the respective quantities. 

The question is whether the linear response follows this 
prediction for the soft modes. The width w of the peak in 
P(a) is, close to the jamming transition, roughly u\\/u± 
because \ctij — tt/2| as uuij/u^^j if uny <C u±jj. It turns 
out that the scaling behavior (fTS")) is consistent with the 
width w of the peak of P(a) for shear deformations, but 
not for compression. There the peak of P{a) does not 
grow as much, and a substantial shoulder for large a re- 
mains even close to jamming: the tendency for particles 
to move towards each other remains much more promi- 
nent under compression. 

Scaling o/wy and u± — The scaling of the distributions 
of U|| and u± has also been probed. The key observation 
is that in Eq. ([8]) the terms ~ U|| and u± have opposite 
signs. What is the relative contribution of these terms, 
and can we ignore the latter? Surprisingly, even though 
S <C 1, Eq. (fTS"|) predicts that the two terms are of equal 
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magnitude in soft modes, and so for linear response one 
needs to be cautious. 

It has become clear that the balance of the terms is 
never so precise as to qualitatively change the magni- 
tude of the energy changes: AE and ^J2ij^ij( u f| %j) 
scale similarly pnl |62| . Hence, the typical values of u\\ 
under a deformation are directly connected to the corre- 
sponding clastic modulus: For compression u\\ is essen- 
tially independent of the distance to jamming (ity ~ e), 
while for shear, u\\ y e Ar/> 1//4 , where e is the magnitude 
of the strain [3l|,[62|. 

The scaling for the amount by which particles in 
contact slide past each other, is more subtle. Numer- 
ically, one observes that for shear deformations, uj_ ~ 
e A(/> -1//4 . The two terms ~ u\\ and ~ wj_ become com- 
parable here, and the amount of sideways sliding under 
a shear deformation diverges near jamming (30l . 13X1 . |62| . 
For compression there is no simple scaling. Combining 
the observed scaling for uy with Eq. (|15[) . one might have 
expected that uj_ ~ e However, the data sug- 

gests a weaker divergence, close to Ar/>~ - 3 . Hence, con- 
sistent with the absence of simple scaling of the peak of 
P(ot) for compression, the two terms oc uii and cx uj_ 
do not balance for compression. Nevertheless, both un- 
der shear and compression, the sliding, sideways motion 
of contacting particles dominates and diverges near jam- 
ming. 



(a) Effective medium theory 



(b) Packings 




K/k 



(c) Random networks 
K/kt 




FIG. 14: Schematic comparison of the variation of shear (G) 
and bulk (K) elastic moduli as function of distance to jam- 
ming, (a) In effective medium theory, all elastic moduli are 
simply of the order of the local spring constant k, and more- 
over, the theory does not account for whether the packing 
is rigid or not. (b) In jammed packings of harmonic parti- 
cles, the bulk modulus K remains constant down to the jam- 
ming transition, where it vanishes discontinuously, whereas 
the shear modulus G vanishes linearly in Az. (c) In random 
networks of elastic springs, both elastic moduli vanish lin- 
early with Az. (From [62( | — Copyright by the Institute of 
Physics) . 



5. Effective Medium Theory, Rigidity Percolation, Random 
Networks and Jammed Systems 

In 1984, Feng and Sen showed that elastic percolation 
is not equivalent to scalar percolation, but forms a new 
universality class [64| . In the simplest realization of rigid- 
ity percolation, bonds of a ordered spring network are 
randomly removed and the elastic response is probed. 
For such systems, both bulk and shear modulus go to 
zero at the elastic percolation threshold [95[ , and at this 
threshold the contact number reaches the isostatic value 
2d [65| . Later it was shown that rigidity percolation is 
singular on ordered lattices (66j . but similar results are 
expected to hold on irregular lattices. 

While it has been suggested that jamming of friction- 
less spheres corresponds to the onset of rigidity percola- 
tion [59], there are significant differences, for example 
that the contact number varies smoothly through the 
rigidity percolation threshold but jumps at the jamming 
transition 0. Nevertheless, it is instructive to compare 
the response of random spring networks of given contact 
number to those of jammed packings — note that the 
linear response of jammed packings of particles with one- 
sided harmonic interactions is exactly equivalent to that 
of networks of appropriately loaded harmonic springs, 
with the nodes of the network given by the particle cen- 
ters and the geometry and forces of the spring network 
determined by the force network of the packing. 

In Fig. (fT4")h a schematic comparison of the variation 



of the elastic moduli with contact number in effective 
medium theory, for jammed packings and for random 
networks is shown. This illustrates that EMT predicts 
that the elastic moduli vary smoothly through the iso- 
static point and that the moduli arc of order of the lo- 
cal spring constant k. This is because effective medium 
theory is essentially "blind" to local packing considera- 
tions and isostaticity. Thus, besides failing to capture 
the vanishing of G near jamming, its prediction for the 
bulk modulus fails spectacularly as well: it predicts finite 
rigidity below isostaticity. Clearly random networks also 
fail to describe jammed systems, as for random networks 
both shear and bulk modulus vanish when z approaches 
Ziso — from the perspective of random networks, it is the 
bulk modulus of jammed systems that behaves anoma- 
lously. 

By comparing the displacement angle distributions 
P(a) of jammed systems and random networks un- 
der both shear and compression, Ellcnbroek et al. con- 
clude that two cases can be distinguished [63]. In the 
"generic" case, all geometrical characterizations exhibit 
simple scaling and the clastic moduli scale as Az — this 
describes shear and bulk deformations of randomly cut 
networks, as well as shear deformations of jammed pack- 
ings. Jammed packings under compression form the "ex- 
ceptional" case: the fact that the compression modulus 
remains of order k near jamming is reflected in the fact 
that various characteristics of the local displacements do 
not exhibit pure scaling. 





FIG. 15: Part of a packing of frictional discs in two dimensions 
for low pressure, zero gravity and friction coefficient fi = 10. 
For this packing, z m 3.06 and 4> ~ 0.77 (this density includes 
rattlers, which are not shown in this image and occur in the 
"holes"). Lines indicate the strength of the normal forces 
— note the large number of near contacts (pairs of particles 
appearing to touch but not connected by a force line). Disc 
color indicates local contact number, clearly identifying the 
large fraction of particles with two contacts only — these do 
not arise in frictionless systems. 



F. Conclusion 

For packings of soft frictionless spheres and in the limit 
of large systems, contact number, packing density, par- 
ticle deformation and (for given force law) pressure are 
all directly linked and at point J the system becomes iso- 
static. The jamming transition for frictionless spheres 
exhibits a number of non-trivial scaling behaviors, all in- 
timately linked to the non-trivial square-root scaling of 
the excess contact number with distance to the isostatic 
jamming point. We have stressed the viewpoint that ge- 
ometry and mechanics are intimately linked for these sys- 
tems, and that near point J, local non-affinity and global 
anomalous mechanical scaling go hand in hand. 



IV. JAMMING OF FRICTIONAL SPHERES 

Here we discuss the rich phenomenology of jamming 
of frictional soft spheres. The crucial difference with 
the frictionless case is that both the packing density 4> c 
and contact number z c at jamming are not unique: both 
depend on the friction coefficient and on the history of 
the packin g, a nd are lower than for frictionless spheres 

[22L m, Mm, Hi M Iz3 — see fi g ure E3 

Jamming and isostaticity no longer go hand in hand 
for frictional spheres. The contact number at jamming, 
z c , can range from d + 1 to 2d, where d + 1 is the iso- 
static value z? for frictional spheres (see section [TV A II 
and Appendix |A"|) . It appears that z c approaches 
only in the limit of fi — > oo and very slow equilibration 
of the packings [69|, [TO, [7l| — see section IIV A 41 In 
all other cases, the number of contacts at jamming is 



FIG. 16: Frictionless (a) and frictional (b) disc in a groove 
[7^ |. (a) In the frictionless case, the system is isostatic, and 
the contact forces (black) balancing the gravitational force 
(blue) are unique, (b) In the frictional case, the system is hy- 
perstatic: contact forces in hard frictional systems are, in gen- 
eral, under-determined. In this example, there are four force 
degrees of freedom (two normal and two frictional forces), and 
only three balance equations (total force in x and y direction 
and torque balance). This leads to a family of solutions (three 
examples indicated in red, orange and green) that balance the 
gravitational force (blue). Which of these is realized depends 
on the history of the system. 



larger than the minimal number needed for force balance 
and rigidity, and frictional packings of soft spheres at 
jamming (or, equivalently, frictional rigid spheres), are 
hyperstatic: z c > z^ so . Hyperstaticity implies that for 
packings of rigid, frictional spheres, the contact forces 
arc not uniquely determined by the packing geometry, as 
was the case for the isostatic packings of rigid, friction- 
less spheres 0, Ef| . An explicit example of this so-called 
indeterminacy of frictional forces is shown in Fig. \TE\ 72% . 

What does the deviation of the critical contact number 
from the isostatic value imply for the scaling of quantities 
such as G, K and uj*1 We will show that these scale with 
distance to the frictional isostatic point, z — z^ so . Thus, 
when the jamming transition is approached, bulk quanti- 
ties in general do not exhibit scaling with distance to the 
jamming point, since, at jamming, z approaches z c > z^ so 
lH, [6?], [H, [7(| IZllj [ill • Scaling with distance to jamming 
can only occur when z c = z? so . Hence, jamming is not 
critical for frictional systems: power law scaling of bulk 
quantities with distance to jamming is the exception, not 
the rule. 

The jamming scenario for frictional soft spheres is de- 
tailed below. We briefly discuss the frictional contact 
laws in section HV II In section HV Al we discuss the prop- 
erties of frictional sphere packings at the jamming thresh- 
old, or equivalently, packings of undcformable frictional 
spheres. We focus on the variation of the range of contact 
numbers and densities as function of /x in sections IIV A It - 
IIV A 31 Finally in section HV A 4l we introduce the concept 
of generalized isostaticity, which is relevant for frictional 
packings that have fully mobilized contacts. Section llVBI 
concerns frictional packings at finite pressures and we 
discuss the (breakdown) of scaling with distance to jam- 
ming. 
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1. Frictional Contact Laws 

Friction is taken into account by extending the con- 
tact force model to account both for normal forces F n 
and tangential forces F t . In the simple Coulomb pic- 
ture of friction, contacts do not slide as long as the ratio 
of tangential and normal forces remains smaller than or 
equal to the friction coefficient /i: |F t |/F n < /i. which 
introduces a very sharp nonlinearity in the contact laws. 
Typical values for fj, relevant in experiments range from 
0.1 to 1, which is where properties of frictional packs vary 
strongly with \i. 

Frictional forces do not only depend on the relative 
position of the contacting particles, but also on their 
history [H H3, IzH [Zi [Za Izi- This is encoded in 
the widely used Hertz-Mindlin model for frictional three- 
dimensional spheres, which takes the normal force F n ~ 
<5 3 / 2 with S the overlap between particles, while the tan- 
gential force increment dF t ~ 5 1 ' 2 dt where dt is the rela- 
tive tangential displacement change, provided F t < /j,F n 
[22l . Fn] . ItH . Studies of friction can also be performed for 
other contact laws, most notably, the linear model for 
which F n ~ S, so that the stiffness of the contacts in the 
normal and tangential direction are independent of the 
normal force and do not vary with distance to jamming 

EH. 
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FIG. 17: (a) Example of the variation of the zero pressure 
contact number z c in two-dimensional rigid discs as function 
of n, smoothly interpolating between the isostatic limits 2d 
(red) for zero friction and d + 1 (blue) for frictional contacts. 
The arched area indicate combinations of contact numbers 
and \x that, while they are not reached in these numerics, are 
perfectly possible — see text (adapted from [68| — Copy- 
right by the American Physical Society), (b) State diagram 
for frictional spheres. While the Random Close Packed, iso- 
static packings obtained for zero friction are compatible with 
all values of /1, a range of packings with lower densities and 
contact numbers open up when /i > 0. For a given prepara- 
tion protocol, there might be a well-defined density (dashed 
curve). Whether there is a well defined lowest packing frac- 
tion for given p, which would define Random Loose Packing, 
is an open question, and the question what the contact num- 
ber of such states would be is open as well ( adapted from 

M). 



A. Frictional Packings at Zero Pressure 

1. Contact Number 

How can the counting arguments for the contact 
number at zero pressure be extended to the frictional 
case? On the one hand, the requirement that contacting 
spheres precisely touch is the same as for the frictionless 
case, and gives zN/2 constraints on the dN particle co- 
ordinates, leading to z < 2d. On the other hand, for fric- 
tional packings, the constraint counting for the zdN/2 
contact force components constrained by dN force and 
d{d — l)N/2 torque balance equations (see Appendix |A")) 
gives z > d + 1, where the isostatic value z£ = d + 1. 
Combining these two bounds, frictional spheres can at- 
tain a range of contact numbers: d + 1 < z c < 2d (see 
Appendix [A"|) . 

It is important to stress that neither bound is sensi- 
tive to the value of /i. What mechanism (if any) selects 
the contact number z c of a frictional packing at jamming? 
The first additional ingredient to consider is the Coulomb 
criterion that for all contact forces \F t \/F n < fi. So, while 
constraint counting allows force configurations that sat- 
isfy force and torque balance for z arbitrarily close to 
z^ so , such configurations arc not guaranteed to be com- 
patible with the Coulomb criterion, and in particular for 
small [i they generally will not be. This is consistent with 
the intuition that a small increase of /1 away from zero 
is not expected to make z c jump from 2d to d + 1. In 
section HV A 41 we will discuss an additional bound on z 



as function of /1. 

Simulations show that in practice z c is a decreasing 
function of ^, approaching 2d at small /i and approach- 
ing d + 1 for large friction coefficient [ig, [68|, IrM l7dl. l7l| 
(Fig. [TTa). However, z c (fi) cannot be a sharply defined 
curve unless additional information about the prepara- 
tion history is given: From the non-sliding condition 
\Ft\/F n < 11 it follows that a packing which is stable 
for a certain value of /u remains so for all larger values 
of /i — increasing the friction coefficient only expands 
the range of allowed force configurations (and does not 
change any of the contact forces). Hence, a numerically 
obtained curve z c {n) at best is a bound for the allowed 
combinations of z c and /1 (see Fig. [TTa). History is a 
second additional ingredient to consider [33], although 
it is remarkable that several different equilibration algo- 
rithms appear to giv e very similar estimates for z c (p) 



2. Density 

The existence of a wide range of statistically different 
frictional packings is also reflected in packing densities, 
which experimentally are more easily observed than the 
contact number. It is well known that packings of spher- 
ical hard particles under gravity (in other words, fric- 
tional spheres close to jamming) can be compacted over 
a range of densities [76|. Different packing densities of 
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these systems do not correspond to deformations of the 
particles, but to changes in the organization of the par- 
ticles. Hence, at jamming, the range of packing densities 
does not go to zero for frictional particles. 

The relation between density and friction coefficient 
can be summarized in a simple state diagram (Fig. 117b ). 
which stresses that random close packing (RCP) is in- 
dependent of /Lt, while the random loose packing (RLP) 
density depends strongly on /it, thus connecting random 
close packing, random loose packing and value of the fric- 
tion coefficient [H, M, M, M, III HI- This dia- 
gram further suggests that the packing density at point 
J may also be seen as random loose packing of friction- 
less spheres (since for = 0, one expects RCP and RLP 
to coincide) — it is the loosest possible packings, rather 
than the densest possible ones, that arise near jamming. 
It should be noted that the definition of RLP is even 
more contentious than RCP, and the debate is wide open 

0,Blzi[zl- 



3. Scaling with fi 

One may now also wonder how the contact number and 
packing density at jamming scale with fx. Qualitative ev- 
idence for scaling was found by Silbert et al. in numerical 
studies of frictional packings (Fig. 2 and 3 from [ill). By 
focussing explicitly on a single preparation protocol, such 
as slow equilibration, this becomes a well posed ques- 
tion loading to the concept of generalized isostaticity, 
defined below. Data for generalized isostatic packings 
suggests that both contact number and density exhibit 
power law scaling with (i for small friction, while for large 
friction, excess contact number and density (defined with 
respect to the infinite friction limit) are also related by 
scaling, although clearly more work is needed to establish 
these scalings firmly [70, [z|[ . 

4- Generalized Isostaticity 

Here we will discuss the role of the frictional forces in 
some more detail, and in particular focus on frictional 
packings for which a large number of contacts are fully 
mobilized, meaning that the frictional forces are maxi- 
mal: \F t \/F n = /i). These packings arise in numerical 
studies when packings are equilibrated slowly for a wide 
range of values of fx. 

The mobilization, m, of a contact is defined as the ra- 
tio \F t \/ (fiF n ), and ranges from zero to one (fully mobi- 
lized) . Earlier numerical data suggested that m generally 
stays away from 1, and that in the limit of large /i, the 
distribution of the mobilization P(m) becomes indepen- 
dent of li [13, EH- Later it became clear that P(m) can 
depend strongly on the preparation history [69j| . Futher- 
more, frictional two-dimensional packings which are very 
slowly equilibrated yield packings for which a substantial 
amount of the contact forces are fully mobilized, meaning 
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FIG. 18: Generalized isostaticity plot, comparing the frac- 
tion of fully mobilized contacts per particle, n m , to the con- 
tact number, z. Data points (open symbols) are for two- 
dimensional systems and for fi ranging from 0.001 to 1000 at 
finite P. The black squares are the corresponding n m and 
z extrapolated to P — 0. The left and bottom axes refer 
to the numerical values for contact number and number of 
fully mobilized contacts per particle, n m , for this specific two- 
dimensional example, while right and top axes give the corre- 
sponding general expressions for higher dimensions. The red 
line denotes the generalized isostaticity line where the number 
of fully mobilized contacts is maximized: n m = d(z — d— l)/2. 
The area to the right of this line refers to generalized hyper- 
static packings, while the area to the left of the red line is 
forbidden (adapted from [7(| — Copyright by the American 
Physical Society). 

that \F t \/F n = n [H,[73,[8l|. One imagines that during 
equilibration, many contacts slowly slide, and when the 
packing jams many contacts are still close to failure — 
such packings are marginal with respect to lowering ll. 

For packings with fully mobilized contacts, the count- 
ing arguments need to be augmented, since at fully mo- 
bilized contacts, the frictional and normal forces are no 
longer independent [70| . Defining the number of fully 
mobilized contacts per particle as n m , the constraints for 
the zdN/2 force degrees of freedom then arc: dN force 
balance equations, d(d—l)N/2 torque balance equations, 
and n m N constraints for the fully mobilized contacts. 
This yields the following relation between z, z? a = d + 1 
and n m [96|: 

z ~ z? so > 2n m /d . (16) 

Surprisingly, for sufficiently slowly equilibrated pack- 
ings and for all values of fi, the values for n m and z 
tend to satisfy this bound when P is lowered to zero 
(Fig. I18|) . Such packings which maximize their number 
of fully mobilized contacts have been referred to as "gen- 
eralized isostatic" packings [zfj, [8l|. These should be 
widely occurring, since most preparation algorithm tend 
to equilibrate slowly 0, HI M, \M, HI [zO, IzH • 

For fully mobilized packings, the amount of fully mo- 
bilized contacts vanishes in the limit of infinite friction 
(see Fig. IT5|) . consistent with the observation that there 
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FIG. 19: Scaling of contact number, w* and elastic mod- 
uli for frictional discs, interacting through three-dimensional 
Hertzian-Mindlin forces, (a) The zero pressure contact num- 
ber, zj does not reach the isostatic limit (z — 3) unless (i is 
very large, (b) The excess coordination number z — zj scales 
linearly with P 1 ' 3 ~ \/A<j!>. (c-d) The characteristic frequency 
of the DOS, to* , scales similarly to z — 3 (e) The bulk 

modulus K (red curves) approaches a plateau for small P, 
while G appears to scale as z — 3 [9(j. (f) As in frictionless 
spheres, the ratio G/K scales with distance to the isostatic 
point, now given by z — z? so — z — 3 (adapted from [Tl|] - 
Copyright by the American Physical Society ). 



z « d + 1. The number of fully mobilized contacts 
is maximal for vanishingly small friction (which we re- 
fer to as yU = + ), where, by continuity, z w 2d, and 
n m ~ d(d — l)/2. Taking into account that each contact 
is shared by two particles, the fraction of fully mobilized 
contacts is (d — l)/2 — hence in two dimensions, 50% 
of all contacts are fully mobilized, in three dimensions, 
100% of the contacts would be fully mobilized for fi = + , 
and in higher dimensions one cannot reach generalized 
isostaticity for /i = + . 

By itself, the inequality (fTB"|) is not a stricter bound on 
z than the ordinary condition z > d+1, since n m (fj,) is un- 
known. However, if we could determine n m (/x), we would 
immediately obtain the bound z = d + 1 + 2 n m (p) / d. It 
is, at present, an open question how n m (fi) can be esti- 
mated or obtained numerically other than through direct 
numerical simulations. 



B. Frictional Packings at Finite Pressure 

Once a mechanically stable frictional packing has been 
created, its linear mechanical response is given by the 
dynamical matrix. For Hertz-Mindlin type interactions, 
each contact can be thought of as being given by two 
springs (one parallel to the contact vector ry, one per- 
pendicular to the contact vector) , the spring constants of 
which are set by the normal force and the Poisson ratio 

Various authors have found that, for essentially all val- 
ues of /i, the excess contact number z — Z c g rows as 
a square root with the excess density [67], El, [Z^ " 
for Hertzian contacts, this is equivalent to stating that 
z — z c ~ P 1 / 3 . However, z c differs from the frictional 
isostatic value d + 1, so that z — z[ l BO does not scale with 
pressure (see Fig.rr5k.bh Note that the slope in Fig.[T5b. 
which represents the prefactor zq in a scaling law of the 
form z — z c — zoV <t> — 4>c does not appear to vary strongly 
with ii. As is the case for frictionless particles, it is es- 
sential to remove rattlers for the count of the contact 
number, but include them for the estimate of the density 
to obtain the square root scaling of z — z c over an appre- 
ciable range [70j. This square root scaling is intriguing 
and, as far as we are aware, without explanation. 

The deviations of z c from the isostatic value imply 
that packings near the (un)jamming transition do not ap- 
proach isostatic packings, and consistent with this, there 
is, in general, no scaling of the mechanical properties as 
function of distance to jamming. 

The mechanical properties do, however, scale with the 
distance to the isostatic point, as measured by the con- 
tact number. First, calculations of the characteristic 
frequency uj* from the density of vibrational states for 
two-dimensional frictional packings show that the vari- 
ation of uj* with /j, and distance to jamming (as mea- 
sured by the pressure P) is very similar to that of z 
(Fig. rr5b). In fact, when this data is replotted as func- 
tion of z — Zjg = z — 3 one finds a linear relation be- 
tween uj* and z — z^ so (Fig.fTTilO. Second, the ratio of the 
shear and bulk modulus exhibits the same phenomenol- 
ogy: G/K scales linearly with z — z? BO = z — 3 (fT9e -f) 
[zlll^|. These findings suggest that, in general, scaling 
is governed by the distance to isostaticity, rather than 
the distance to jamming. 

The contact number and geometry of the packings 
change smoothly with /i ffH , [73[ , while the mechanical 
behavior exhibits a discontinuous jump from n = to 
/j = + . This is caused by the fact that when fric- 
tion is included, the nature of the dynamical matrix 
changes completely, because the tangential contact stiff- 
nesses jump from zero to a finite value. When the tangen- 
tial stiffness is varied smoothly from \i = to finite fric- 
tion, the mechanical properties vary smoothly also [53[ . 

Finally a word of caution regarding the notion of gener- 
alized isostaticity and the role of fully mobilized contacts 
for scaling away from jamming. In the calculations pre- 
sented above fully mobilized contacts are treated as or- 
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dinary elastic contacts. Strictly speaking, such marginal 
contacts cause a breakdown of linear response. One may 
argue that tiny perturbations would simply let the fully 
mobilized contacts relax to almost-fully-mobilizcd, after 
which linear response would no longer be problematic. 
Taking the opposite view, Henkes et al. have recently 
shown that if the dynamical matrix is calculated under 
the assumption that fully mobilized contacts can slide 
freely, the characteristic frequency u>* scales and vanishes 
with P for all values of n — provided one considers sys- 
tems that approach generalized isostaticity [53[ . 



C. Conclusion 

Jamming of frictional grains can be seen as a two-step 
process. The first step is the selection of a contact num- 
ber, z, given the friction coefficient, pressure and proce- 
dure. In the second step, in which the mechanical proper- 
ties of the packing are determined, everything scales with 
z — z^ so . The crucial difference with frictionlcss spheres 
is that the contact number z c at the P = jamming 
point in general does not coincide with z- 1 . Most quan- 
tities are governed by the contact number and scale with 
distance to isostaticity, while the contact number itself 
scales with distance to jamming. 



V. JAMMING OF NON-SPHERICAL 
PARTICLES 

New phenomena occur in packings of non-spherical 
particles, and here we briefly discuss the jamming sce- 
nario for frictionlcss ellipsoids. 

First, configurations for hard (or zero pressure) fric- 
tionlcss ellipsoids pack more densely and have larger con- 
tact numbers than frictionless spheres [H, [49|, y5(| [H[. 
As we discuss in section IV Al both the increase in den- 
sity and in contact number away from the sphere limit 
are continuous but not smooth — plots of cf> and z as 
function of the cllipticity show a cusp at the sphere limit 
(Fig [20]). 

Second, the counting arguments for general ellipsoids 
suggest that at jamming, ellipsoids attain z — z; so = 
d(d + 1). However, weakly aspherical ellipsoids actually 
attain a contact number arbitrarily close to the sphere 
limit z = 2d. As a consequence, (weakly) ellipsoidal 
packings are strongly hypostatic (underconstrained) near 
jamming. This leads to questions about the relation be- 
tween contact number, rigidity and floppy modes (section 

EB). 

Third, the question arises whether quantities such as z 
and to* exhibit scaling, either as function of the pressure, 
as function of the asphericity or as function of distance 
to either the spherical or the ellipsoidal isostatic point — 
the partial answers to these questions, based on recent 
studies of the density of states [ill [5(| will be addressed 
in section IV CI 
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FIG. 20: (a) Packing fraction <f> of spheroids (open symbols) 
and general ellipsoids (blue symbols) as function of the as- 
phericity q. The density shows a cusp at a = 1 (sphere 
limit). The orange line indicates the HCP packing density 
« 0.74, which is almost reached by random packings of some 
ellipsoids, (b) The contact number, z, for the same spheroids 
and ellipsoids as shown in panel (a) also shows a cusp at 
a = 1. The red, green and blue lines at z = 6, z — 10 and 
z = 12 indicate the isostatic contact numbers for spheres, 
spheroids and ellipsoids (adapted from [83l ] — Copyright by 
The American Association for the Advancement of Science). 



Packings of Spherocylinders, Spheroids and 
Ellipsoids 



Spherocylinders — Early indications for surprisingly 
dense packings of non-spherical particles come from stud- 
ies of spherocylinders, particles consisting of a cylinder of 
length a and diameter 1, which on both ends are capped 
by a half sphere. For zero a, these are spheres, while 
the large a limit is relevant for the loose packings of thin 
(colloidal) rods 84]. Williams et al. studied the packing 
fraction and contact numbers of such spherocylinders nu- 
merically, and found that both the packing fraction <j> and 
contact number z increase when a is increased, reach a 
maximum for a ~ 0.4, and then decrease [85]. The den- 
sity peaks at a value of 0.695, substantially larger than 
the typical values for random close packing of spheres 
~ 0.64, while for large a > 10 the density decays as 
4> ~ 1/a, consistent with arguments given before |84| . 

The contact number in these simulations was found to 
start out at z ss 5.8 for a = 0, and increased until it 
reached z w 9 for a m 0.4. The initial value is close to 
the isostatic number for spheres (6), while the peak value 
is similar to the isostatic number for rods (10) [86j. 
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Spheroids and Ellipsoids — In seminal work, Donev 
et al. explored the packing properties of hard spheroids 
and ellipsoids |83| . As shown in Fig [20k . the density 
of spheroids (axis: 1 : 1 : a) exhibits a cusp-like local 
minimum for the pure spherical case a = 1, and reaches 
two local maxima: oblate (disc-like) spheroids at a w 0.6 
pack at a density <fi w 0.70 and prolate (cigar-shaped) 
ellipsoids at a = 1.5 pack even denser at <j) ~ 0.715. Note 
that the spheroid packing density only drops below the 
random close packing value for spheres for very strongly 
oblate (a < 0.25) or prolate (a > 4) particles. 

Even larger packing densities can be obtained for tri- 
axial ellipsoids, and for the case that the axes are given 
as 1/a : 1 : a, the maximum packing density peaks at 
0.735 for a w 1.25, [H,[83j|. This density is surprisingly 
close to the density ss 0.74 obtained for fee and hep pack- 
ings, which arc the densest possible packings for spheres 
- but those are crystals, whereas the ellipsoidal pack- 
ings do not show any appreciable oricntational ordering. 
Finally, crystals of ellipsoids can be packed even denser, 
with the highest density currently known, 0.7707, is ob- 
tained in non lattice periodic packings of spheroids with 
either a > v3 or a < l/Vo [871 ] . 

The contact number grows monotonically with as- 
phericity, from a value m 2d for the spherical case to 
values close to the corresponding higher isostatic number 
for ellipsoids: The contact number for the spheroids mea- 
sured for strongly oblate or prolate appears to level off at 
values around 9.8 (the corresponding isostatic number is 
10), and for ellipsoids one reaches 11.4 (the corresponding 
isostatic number is 12) [83|. (The contact numbers in the 
disordered ellipsoidal systems are difficult to obtain accu- 
rately from numerics, in particular for hard particles — 
since, similar to hard spheres, one expects anomalously 
many near contacts [42|, [47j]). It is noteworthy that the 
contact numbers reach these asymptotic values at the 
same asphericities where the density is maximal. Recent 
work on two-dimensional ellipses 0, |49[ and three di- 
mensional spheroids [5(| confirm these trends in contact 
numbers. 



B. Counting arguments, Floppy modes and 
Rigidity of Ellipsoids 

The counting arguments for general ellipsoids suggest 
that at jamming, ellipsoids attain z = z- lso = d(d + 1). 
However, weakly aspherical ellipsoids actually attain a 
contact number arbitrarily close to the sphere limit z = 
2c?. Hence counting arguments suggest that packings 
of weakly ellipsoidal particles possess a large number of 
floppy modes. Are these packings stable? 

As a first step in understanding such packings, it is 
helpful to think about weakly aspherical ellipsoids that 
approach the sphere limit. The number of floppy modes 
in such an undcrconstraincd system equals (A r /2)(zj SO — 
z), which for the sphere limit (where z — > 2d) equals 
N(d(d— l)/2). What are these floppy modes? 



The key observation is that in the counting arguments 
for ellipsoids, the rotational degrees of freedom are taken 
into account while for spheres, where they correspond 
to trivial rotations of the particles, these are ignored. 
When these rotational degrees of freedom are also taken 
into account for frictionless spheres, one obtains precisely 
N(d(d— l)/2) trivial floppy modes, corresponding to the 
trivial rotational deg rees of freedom of individual friction- 
less spheres [HI, l50l | These floppy modes do not affect 
the rigidity of the packings, which suggests that, in gen- 
eral, absence of floppy modes may be a sufficient but not 
a necessary condition for rigidity [48| . 

From the perspective of constraint counting of the 
contact forces, something similar happens in the sphere 
limit: how do dN force degrees of freedom satisfy both 
dN force balance equations and also all the additional 
torque balance equations? The answer is simple: for fric- 
tionless spheres, the torques exerted by each contact force 
is zero, and so torque balance is trivially satisfied. 

The key question, however, is what happens to hypo- 
static packings at finite asphericity and pressure. The 
full answers are not known, but two recent studies on 
the density of vibrational states for soft frictionless bidis- 
perse two-dimensional ellipses j49| and three-dimensional 
spheroids [5(| provide important ingredients that we will 
discuss below. 



C. Jamming of Ellipsoids 

The main findings for the density of states of ellip- 
soidal particles are shown in Fig. [5T] Close to the sphere 
limit, where the contact number is far below the rele- 
vant ellipsoidal isostatic value, the density of states con- 
sists of three bands: first, a number of zero frequency, 
floppy modes corresponding to the degree of hypostatic- 
ity, second, a band of rotational modes, and third, a 
band of translational modes, corresponding to the trans- 
lational modes present for the pure sphere case. When, 
for increasing pressure and/or strong ellipticity, the con- 
tact number starts to approach the ellipsoidal isostatic 
value, the rotational and translational bands hybridize 
and merge. Finally, when the contact number exceeds 
the ellipsoidal isostatic value, the floppy modes have van- 
ished and the characteristic frequency of the remaining 
single band density of states scales with distance to the 
ellipsoidal isostatic value. 

The counting arguments provide a clear picture of the 
number of modes per band, as shown in Fig. 1221 where 
the variation of these numbers with contact number is 
shown for the case of spheroids in 3d. 

First, the vibrational modes present for the spherical 
case are only weakly perturbed by the inclusion of weak 
ellipticity, so their number still equals dN. The particle 
motions of modes in this band are essentially transla- 
tional, and the characteristic frequency of this band, u>*, 
still scales with z~z^ clc , noi with z — z ; e so ips . Hence, this 
part of the density of states is smoothly perturbed when 
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FIG. 21: Schematic scenario for the density of states for Mo- 
tionless soft ellipsoidal particles, based on [i^, (a-d) Den- 
sity of states as function of distance to the spherical limit. 
The grey, blue, red and green colors refer to floppy modes, 
rotational modes, translational modes and hybridized modes 
respectively, (a) For frictionless spheres, one usually only con- 
siders the translational band (red), but when one takes the 
rotational degrees of freedom into account, a large number of 
trivial floppy modes occur (dashed grey line), (b) For con- 
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sphere 



the density of states ex- 



hibits three bands, and the characteristic frequencies u a 



scale with asphericity and z — z 



sphere 



respectively 



and 
see 



text, (c) For contact numbers approaching z — z^ p , the rota- 
tional and translational band merge, (d) For contact numbers 
above z — zj ip , there are no floppy modes and the character- 
istic frequency scales with z — zL v 



going from the sphere to the weakly ellipsoidal case. 

Second, for z < z? ip the system is underconstrained, 
and the crucial observation is that here there are [z — 

z iso P )/2 floppy modes. In the sphere limit, these modes 
are the trivial local rotations, and away from the sphere 
limit most of these modes survive and become delocalized 
— their precise nature is not fully understood yet. 

Thirds, at finite pressures and/or finite asphericities, 
(z — z^ ere )/2 modes emerge from the zero frequency 
band and attain finite frequencies. This is the rotational 
band: particle motions of modes in this band are essen- 
tially rotational, and the vibration frequencies are below 
those of the translational band. This allows the defini- 
tion of a characteristic maximal frequency of the rota- 
tional band uj s , which is found to scale with the degree 
of asphericity |1 — a\, but is essentially insensitive to the 
pressure. 

Fourth, for large pressure and asphericity, the contact 
number approaches the relevant ellipsoidal isostatic num- 
ber, the rotational and translational bands start to ap- 
proach each other ( lu* — uj s <C 1), the modes hybridize 
and these two bands eventually merge. In the regime 
where the contact number exceeds the relevant ellipsoidal 
contact number, there are no more floppy modes. The 
only band of vibrational modes then has a mixed trans- 
lational/rotational character, and its characteristic fre- 
quency, uj^ , scales with distance to the relevant ellipsoidal 
isostatic point: 10 s ~ z — zf^ p . 

Finally, note that for weakly elliptical systems that are 
hypostatic, the counting argument implies that the forces 
must be non-generic — one still has more equations of 
force and torque balance than one has force degrees of 
freedom. In terms of the elastic energy landscape, one 
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FIG. 22: Schematic representation of the number of modes 
per band for the specific case of spheroids in three dimensions 
— from [Htj — Copyright by the Institute of Physics 



imagines that near such systems there must exhibit many 
directions in phase space where the second derivative is 
zero (leading to quartic modes (49[), but a deep under- 
standing is missing. 



D. Conclusion 

Jamming of frictionless ellipsoidal particles is surpris- 
ingly similar to that of frictionless spheres, despite the 
strongly hypostatic nature of weakly aspherical packings. 
The crucial observation is that frictionless spheres can 
also be seen as strongly hypostatic near jamming, as they 
possess a large number of trivial floppy modes. Most of 
these modes remain at zero frequency for weakly ellip- 
soidal particles, even though their spatial structure is no 
longer trivial, and these modes do not appear to affect 
the rigidity of packings of frictionless ellipses. 



VI. SUMMARY, OPEN QUESTIONS AND 
OUTLOOK 



The jamming scenario for disordered packings of soft, 
purely repulsive particles at zero temperature and shear, 
as described above, can be seen as a two step process. 
First, for a given pressure, contact law and preparation 
protocol, a packing with a certain contact number, z, 
is created. Second, the mechanical characteristics such 
as elastic moduli and density of states depend on the 
difference between the actual contact number and the 
relevant isostatic value. 

Depending on the particles' friction or shape, the con- 
tact number may span a range of values — see Fig. l23l for 
this range for P — > 0. For frictionless particles it appears 
that the contact number at jamming is independent of 
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FIG. 23: Conjectured range of selected contact numbers at 
jamming, i.e., at P = 0, as function of the friction coefficient 
fi and ellipticity e. The red dot indicates the isostatic limit for 
frictionless spheres at (fj, = 0, e — 0), the green line indicates 
the isostatic limit for frictionless ellipsoids ((i = 0,e / 0) 
and the blue plane indicates the isostatic limit for frictional 
particles (fi 7^ 0). The contact number is precisely selected 
in the frictionless plane, and for sufficiently large ellipticity 
the contact number crosses the isostatic value (green dot). 
Once friction comes into play, a range of contact numbers are 
allowed. For a given fi and e, the upper bound is given by the 
selected contact number for frictionless ellipsoids at /i = 0, 
while the lower bound is given by the generalized isostaticity 
limit — i.e., for finite fi, the maximal number of contacts is 
fully mobilized here, and only for fj, — > 00 does z reach the 
frictional isostatic value z = d + 1. 



the preparation procedure, even for finite pressures. For 
frictional particles, a range of contact numbers arises and 
the history becomes crucial. 

Jamming of frictionless soft spheres constitutes a spe- 
cial case, since here the isostatic contact number (exclud- 
ing the trivial rotational degrees of freedom of the parti- 
cles) is reached at the jamming threshold. The counting 
for ellipsoidal particles takes these rotational degrees into 
account, which leads to strongly hypostatic packings near 
jamming — however, the associated zero modes do not 
appear to contribute to the mechanical properties of the 
packings. Furthermore, the perturbation from spheres 
to weak ellipsoids is smooth, when the trivial rotational 
modes for the spheres is included. 

Friction, however, acts differently. Given a certain 
preparation procedure, the change in contact number is 
smooth with fj,. However, the frictional interactions are 
such that at the level of the dynamical matrix, the in- 
clusion of arbitrary small friction introduces a discontin- 
uous change. For any value of the friction the tangential 



stiffness takes on a finite value which leads to contribu- 
tions to the dynamical matrix of order one, contributions 
which are absent in the frictionless case. Friction become 
a smooth perturbation only when the tangential stiffness 
is varied smoothly with /i. 



A. Open Questions 

A crucial question is that of experimental relevance. 
Many recent predictions of the theory should be observ- 
able in experiment, in particular for frictionless systems 
such as foams and emulsions, but very few have been ob- 
served so far. Frictional packings have been explored the- 
oretically far less than frictionless systems, despite their 
obvious experimental relevance (T3 , |88| . How many dif- 
ferent order parameters does one need to characterize the 
statistics of generic frictional packings? 

More work is needed to clarify the notion of random 
loose packing [z3, [zl, [z!] , and to unravel the role of pack- 
ing protocols. What is the underlying distribution of pos- 
sible contact numbers and densities for frictional spheres, 
given a certain pressure and friction coefficient? Do RCP 
and RLP correspond to sharp gradients in this distribu- 
tion? Are the RCP and RLP limits identical for fric- 
tionless packings? Random packings of spheres are much 
looser then random packings of non spherical particles — 
can we understand why? 

It is, in many cases, unknown how results obtained 
for frictionless spheres extend to more complex systems. 
For example, do a diverging length scale and a singularly 
non-affine response arise when frictional spheres or el- 
lipses approach their isostatic limit(s)? What about the 
elastic moduli [7l|, [HJ? Similarly, what is the jamming 
scenario for more general particles, such as frictional el- 
lipses and non convex particles that may share multiple 
contacts? What is the scenario for more general interac- 
tions (attraction, long range...)? 

Given the central role of the square-root scaling of the 
contact number with distance to jamming, it would be 
useful to probe the connection to the square-root singu- 
larity of g(r) — the argument outlined in section fill D II 
assumes displacements to be primarily affinc, while near 
J, the displacements are singularly non-affine and di- 
verge. What may happen is that the relative displace- 
ment of particles that are not in contact are not strongly 
non-affine — we don't know. For frictional spheres it is 
not understood whether z — z c exhibits true square root 
scaling with excess density, and whether g(r) exhibits 
similar scaling behavior there. 

Essentially all the work discussed above focuses on av- 
eraged quantities and linear response. For finite systems, 
contact numbers, moduli etc exhibit significant differ- 
ences in different realizations 0, [52j ■ Can we understand 
these fluctuations near jamming? What is the nonlinear 
yielding behavior of systems near jamming 0]? 

A whole host of new phenomena arise when jammed 
systems are put under shear stress, and possibly are made 
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to flow d, @, H, Q , or when systems of finite temperature 
[Tol . [IH are considered. Can these phenomena be con- 
nected in a meaningful manner to the zero shear, zero 
temperature limit? 

B. Outlook 

Jamming is cool [l| , as it provides a framework to ap- 
proach the mechanics of disordered systems. The studies 
of the simplest case of static soft frictionlcss spheres have 
demonstrated that such systems exhibit rich spatial orga- 



nization and anomalous mechanical properties near the 
isostatic/jamming limit. Important tasks for the coming 
years include exploring the relevance of these observa- 
tions for experimental observations and for systems with 
more complex interactions. New horizons are emerging 
for systems at finite temperature and in particular for 
flow near jamming — as attested by the rich phenomenol- 
ogy of flowing foams, suspensions and granular media. 
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APPENDIX A: COUNTING ARGUMENTS FOR 
THE CONTACT NUMBER 

By constraint counting one can establish bounds on the 
contact number (43|. First, one may require that floppy 
modes, deformations that in lowest non trivial order do 
not cost energy, are absent. This yields a lower bound 
on the contact number. Packings that violate this second 
constraint are called hypostatic, packings that marginally 
fulfill this constraint are isostatic, and packings that ful- 
fill this constraint are called hyperstatic. 

Note that the same lower bound on the contact number 
is obtained by requiring that all contact forces balance. 
As we will see, this is because the number of independent 
degrees of freedom necessarily to describe changes in the 
energy at a contact equals the number of force degrees 
of freedom per contact. Therefore, the requirement that 
floppy modes are absent is equivalent to the requirement 
that the contact forces balance, and often the counting 
argument that yields the lower bound on z is phrased in 
terms of the contact forces. 

Secondly, for packings at jamming, one arrives at a sec- 
ond constraint, which follows from the requirement that 
the particles are undeformcd at jamming. This yields an 
upper bound on the contact number. Violations of this 
second condition are possible for special (non-generic) 
packings, such as perfect crystals. 

As we will see, for frictionless particles the first and 
second bounds coincide. This does not necessarily im- 
ply that the corresponding contact numbers are real- 
ized at jamming: numerically it is found that frictionless 
spheres are indeed isostatic at jamming [2J, while weakly 
aspherical frictionless ellipsoids are strongly hypostatic 
[IS 5sL Hs| ■ For frictional particles the two bounds never 
coincide, and numerically it is found that frictional par- 
ticles are almost always hyperstatic at jamming. 

Below we present the counting arguments in detail, 
for packings of TV soft particles in d dimensions which 
interact through contact forces, and for which the con- 
tact number z, is defined as the average number of con- 
tacts per particle. Note that the total number of contacts 
equals TVz/2 — each contact is shared by two particles. 
As we will find below, to perform these counting argu- 
ments we need to know the number of force components 
per contact, or cquivalently, the number of independent 
degrees of freedom necessarily to describe changes in the 
energy at a contact (/), the geometrical number of de- 
grees freedom per particle (x) and the number of force 
balance equations per particle (6). 

Absence of Floppy Modes — The counting that 
follows from requiring that there are no floppy modes can 
most easily be carried out by considering AE, the change 
in elastic energy as function of deformation of a certain 
packing. The number of terms contributing to AE equals 
the number of contacts, Nz/2, multiplied with /, the 
number of independent degrees of freedom necessarily to 
describe changes in the energy at a contact. AE is a 
function of all TVc? positional degrees of freedom, and all 



additional orientational degrees of freedom which are not 
symmetries — zero for spheres, 2TV for spheroids in three 
dimensions, and d(d — l)TV/2 for general ellipsoids. We 
denote these number of degrees of freedom relevant for 
AE by b. 

Absence of generic floppy modes requires that the num- 
ber of terms contributing to AE exceeds the number of 
degrees of freedom: zf jl > b. 

For frictionless particles, / equals one, because energy 
changes result from (de)compression of contacts only, 
while for frictional particles, / equals d, since relative 
motion of contacting particles in all directions are rele- 
vant. 

The situation for b is simple for frictional particles, 
where all positional and orientational degrees of freedom 
are relevant and b = d(d+ 1 )/2. For frictionless particles, 
b depends on the symmetries. For frictionless spheres, 
only translational degrees of freedom are important and 
b = d. For frictionless spheroids in three dimensions, two 
additional rotational degrees of freedom come into play 
and 6 = 5, while for general frictionless ellipsoids, all 
rotational degrees are relevant and b = d(d + l)/2. 

Equivalence of Floppy Mode and Force Balance 
Counting — The requirement zf /2 > & is exactly the 
same as requiring that there are sufficient contact forces 
in the system so that they generically can be expected to 
balance: the number of contact force degrees of freedom 
per particles is zf/2 and the number of equations that 
need to be satisfied equals b. The number of relevant 
particle degrees of freedom in the energy expansion thus 
corresponds to the number of force balance equations, 
and the number of terms in AE (= number of constraints 
needed to generically avoid Soppiness) correspond to the 
number of force degrees of freedom — changes in energy 
and forces are directly linked. 

Note that even though the role of constraints and de- 
grees of freedom interchanges when altering the picture 
between absence of floppy modes and satisfaction of force 
balance, so does the requirement (floppy modes: making 
sure there are no generic solutions, force balance: making 
sure there are generic solution), and in the force balance 
picture one ends up with precisely the same inequality: 
zf/2 > b. 

Touch — The conditions that particles precisely 
touch yields Nz/2 constraints on the degrees of freedom 
of the particles. Denoting the number of geometric de- 
grees per particle as x, the condition that for generic 
packings there should be less constraints than degrees of 
freedom yields z/2 < x. 

For the particles that are considered here (spheres and 
ellipsoids with and without friction), the number of de- 
grees of freedom per particles are their d positional co- 
ordinates, to which ellipsoids add their relevant angular 
degrees of freedom. For general ellipsoids, these yield 
d(d + l)/2 degrees of freedom — for spheroids in three 
dimensions (an ellipsoid with two equal axes, which thus 
has one symmetry of rotation — see section [V} these 
yield 5 degrees of freedom. The resulting counting of 
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Particle 


/ 


X 


b 


Touch 


Rigidity 


Range 










z/2 < x 


zf/2 > b 




Frictionless Sphere 


1 


d 


d 


z<2d 


z>2d 


z = 2d 


Frictional Sphere 


(/ 


d 


d(d+l)/2 


z<2d 


z > d+1 


d+1 < z < 2d 


Frictionless Spheroid 


1 


5 


5 


z<10 


z>10 


z = 10 


Frictional Spheroid 


3 


5 


6 


z < 10 


z>4 


4 < z < 10 


Frictionless Ellipsoid 


1 


d(d+l)/2 


d(d+l)/2 


z < d(d+l) 


z > d(d+l) 


z = d(d+l) 


Frictional Ellipsoid 


d 


d(d+l)/2 


d(d+l)/2 


z < d(d+l) 


z > d+1 


d+1 < z < d(d+l) 



TABLE I: Results of "Maxwell" constraint counting for a range of different type of soft particles. As explained in the text, / 
denotes the number of force components per contact, x denotes the geometrical number of freedom per particle and b denotes 
the number of balance equations per particle. 



x and corresponding inequalities are listed in table [U In 
particular, for frictional particles, the lower bound for the 
contact number is d+1, while for frictionless particles it 
depends on the symmetries of the particles. 

Results — The resulting inequalities are listed in ta- 
ble |H Note that the upper bounds for z coincide for fric- 
tional and frictionless particles, as this number only de- 
pends on the geometrical number of degrees of freedom. 



The inequalities can be summarized as follows: For fric- 
tionless particles, / equals one, b = f and the lower and 
upper bounds coincide at z = 2x = 2b/ f. For frictional 
particles, 2i > 2b/ f, the lower and upper bounds do not 
coincide, and a range of contact numbers is allowed at 
jamming. 



